Gate defined zero- and one-dimensional confinement in bilayer graphene 



Augustinus (Stijn) M. Goossens 1 !!! Stefanie CM. Driessen 1 , Tim A. Baart 1 , 
Kenji Watanabe 2 , Takashi Taniguchi 2 , and Lieven M.K. Vandersypen 1 
1 Kavli Institute of Nano science, Delft University of Technology, P.O. Box 5046, 

2600 GA Delft, The Netherlands, 2 Advanced Materials Laboratory, 
National Institute for Materials Science, 1-1 Namiki, Tsukuba, 305-0044, Japan 

(Dated: March 4, 2013) 

We report on the fabrication and measurement of nanoscale devices based on bilayer graphene 
sandwiched between hexagonal boron nitride bottom and top gate dielectrics. The top gates are 
patterned such that constrictions and islands can be electrostatically induced by applying appro- 
priate voltages to the gates. The high quality of the devices becomes apparent from conductance 
quantization in the constrictions at low temperature. The islands exhibit clear Coulomb blockade 
and single-electron transport. 



Confinement of charge carriers in graphene has 
been heavily investigated since _graphene flakes on a 
substrate were first measured. [l|, [2[ Nanopatterning 
graphene into nanoribbons or small islands has been a 
widely used strategy for both zero- and one-dimensional 
confinement. [3KZI Unfortunately, as a result of the etch- 
ing process that is used for nanopatterning, the edges 
of the devices are very rough. This edge disorder, ag- 
gravated by substrate- induced disorder, leads to rather 
irregular device behaviour. For instance no quantized 
conductance was observed in etched constrictions, the 
transparency of barriers varied non-monotonously with 
gate voltage, no atom-like shell- filling has been observed 
in quantum dots and in many cases quantum dots fell 
apart into multiple islands. 

Last year, suspended single layer graphene sheets nar- 
rowed by current induced heating showed quantized steps 
in conductance of 2e 2 /h as the Fermi wavelength was 
varied. [Si Presumably these constrictions were not only 
narrow but also short, as observed in recent transmission 
electron microscope measurements [9], reducing the ef- 
fects of edge disorder. However, the formation of these 
constrictions is hard to control, making it difficult to re- 
produce these results. Clean ribbon edges can be ob- 
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tained reproducibly by unzipping carbon nanotubes, and 
a well-behaved quantum dot formed in such a nanorib- 
bon was recently reported. [To| Barriers were formed at 
metal Schottky contacts, but such barriers are not tun- 
able, limiting follow-up work. Moreover, as the ribbons 
are not obtained lithographically but are dispersed from 
solution, they face many of the limitations of the carbon 
nanotube they originate from. 

The ideal device would confine charge carriers in the 
bulk, far from (disordered) edges, have well controlled 
tunnel barriers, and enjoy all the design freedom offered 
by lithography. All these requirements can be satisfied 
using patterned electrostatic gates, provided a band gap 
is present. This calls for the use of bilayer graphene 
rather than monolayer graphene, as in bilayer s a band 
gap can be induced by an electric field perpendicular 
to the layers. [Tll - fl3( | An additional requirement for clean 
confinement is to minimize substrate-induced disorder. 
Substrate disorder can be eliminated in suspended de- 
vices with suspended top gates, as shown in recent work 
by Allen et al.[14] However, it would be highly desirable 
to realize devices of comparable quality on a substrate, as 
this would facilitate integration of complex devices. Cur- 
rently the cleanest gate dielectric available for graphene 
devices is hexagonal boron nitride (hBN), and mobilities 
reported on such substrates approach those of suspended 
devices. flBj 

Here we report 1-dimensional and 0-dimensional elec- 
trostatic confinement in bilayer graphene sandwiched in 
hBN dielectrics. The devices have split top gates and a 
global back gate, which we bias so that a gap is opened. 
The Fermi level is tuned inside the gap in the regions 
below the top gates so that they become insulating. De- 
pending on the combination of top gates that is biased, 
we form 1-dimensional channels or 0-dimensional islands. 

A schematic and AFM image of a device is shown in 
Fig. [TJ We first deposit a 14 nm thick hBN flake by 
mechanical exfoliation on a silicon wafer coated with a 
silicon oxide (Si02) layer of thickness tsio 2 — 285 nm. 
On top of the hBN we transfer a bilayer graphene flake 
(~ 22/im long and ~ 3/im wide, its bilayer nature was 
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confirmed by Raman spectroscopy) using a dry transfer 
method following the protocol of [15| (at a temperature 
of 100 °C to remove any water absorbed on the surface 
of the graphene and hBN flakes). The sample was subse- 
quently annealed in an oven at 400 and 450 °C (Ar 2400 
seem, H2 700 seem) to remove residues induced by the 
transfer process. Cr(5 nm)/Au(95 nm) electrodes are 
fabricated using electron-beam lithography (EBL). We 
annealed the samples again (same flow rate as the first 
annealing step, T = 300, 350 and 440 °C) to remove fab- 
rication residues. This did not give the desired sample 
quality. Hence we applied the recently developed me- 
chanical cleaning technique [16| followed by dry transfer 
of a 50 nm thick hBN flake, which will act as top gate di- 
electric. In a two-step EBL process we deposited Cr/Au 
top gates. We defined several gate patterns between the 
contacts on this flake and report here on two 'quantum 
dot' top gate structures, one with a lithographic diameter 
of 320 nm (device A) and one of 250 nm (device B). The 
separation between the top gates that together define a 
barrier is less than 30 nm. For device B, TGI and TG2 
were unintentionally connected. 

This specific graphene sample was cooled down multi- 
ple times. The maximum overall mobility measured at 
T=35 mK was ~ 36000 cm 2 /Vs (four-terminal con- 
figuration). The overall mobility measured during the 
last cool-down was much lower ( ~ 6000 cm 2 /Vs, two- 
terminal configuration). This degradation of electronic 
quality was caused by exposure to an electron beam dur- 
ing imaging in a scanning electron microscope. Never- 
theless we believe that the graphene sandwiched in hBN 
remained of high quality. 

We set the back gate to a large negative voltage and 
tune the top gate to a voltage that compensates for the 
doping induced by the back gate. For typical values 
Vbg = — 50 V and Vtg = 9 V the displacement field 
is D ~ 0.6 V/nm, which translates into a theoretically 
predicted band gap of ~ 50 meV. G3, Ell As in earlier 
work, the transport gap is substantially smaller 
but still large enough to realize quantum confinement, as 
we will see. 

In Fig. [2^l we show a top gate trace at a large nega- 
tive back gate voltage, taken at low temperature. Clearly 
visible is a region of suppressed conductance, with a re- 
markably clean transition region on the hole side, and a 
somewhat less clean transition and lower conductance on 
the electron side. This behavior can be expected since the 
leads are p-doped by the backgate, and a pnp-junction is 
formed on the electron side. The pnp-junction decreases 
the transparency of the device as the charge carriers have 
to Zener tunnel through the induced band gap. On the 
hole side the conductance suppression with top gate volt- 
age is very smooth and well behaved. This is in strong 
contrast to similar top gate sweeps for graphene nanorib- 
bons, which exhibit very irregular pinch-off characteris- 
tics when the Fermi energy is swept into the transport 
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Zooming in on the steep flank of the pinch-off curve, 
we observe several plateaus in conductance, with a value 
that is independent of Vbg (Fig. Wp)- The traces have 
been corrected for filter resistance, current amplifier in- 
put resistance and for a background resistance (Rb) con- 
sisting of contact resistance, Maxwell spreading resis- 
tance and graphene lead resistance. The last two terms 
are dependent on back gate voltage. For each value of 
Vbg we subtracted a value of Rb comparable to the two- 
terminal resistance (ifer) of the graphene flake with all 
top gates at ground (see inset fig. [2)3), fine tuning Rb 
such that the conductance G at the upper plateau is 
6e 2 //i.0 The other two plateaus consequently appear 
at G = 4e 2 /h and G = 2e 2 /h. The same sequence of 
steps was observed for device A. (in device B, a less well 
developed feature can be seen just below G = 3e 2 //i, 
which does not appear in device A; its origin is unclear) 
This behavior is consistent with transport through one- 
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FIG. 2. (a) Conductance versus the voltage on TGI and TG2 
(device B, T = 35 mK, V bias = 0.5 mV, V bg = -49.74 V). The 
other top gates of the device are set to 5 V, far from pinch- 
off, so that we are dealing with only one constriction. All 
measurements presented in this paper are conducted in two- 
terminal DC voltage bias configuration. The conductance has 
been corrected for filter resistance and current amplifier input 
resistance. In the inset we show cartoons of the bandstructure 
for pp'p and pnp measurement regimes. Solid lines indicate 
valence and conduction band edges and dashed line the po- 
sition of the Fermi level, (b) Top gate traces for device B 
at back gate voltages from -49.26 V (leftmost) to -49.74 V 
(rightmost) in 10 equal steps (Vbias = 0.5 mV, T = 35 mK). 
The solid green line in the inset is the subtracted Rb- The 
dashed blue line is the two-terminal resistance of the graphene 
flake with all top gates at ground, (c) Channel width W n for 
modes n = 1, 2 and 3 for device A (blue triangles) and for two 
measurements on device B (green circles and squares). The 
widths are calculated in the square- well potential approxima- 
tion W = f A/. The Fermi wavelength at the n-th plateau 



is estimated via A/ 



Ctg is the top gate 



capacitance per unit area, extracted from the slope of the 
position of the v — 4 plateau in top gate voltage against mag- 
netic field (Fig. EJd). AVtg,™ is the difference between the 
top gate voltage in the middle of the n-th plateau and the 
top gate voltage at charge neutrality. 
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FIG. 3. (a) Conductance versus TGI and TG2 voltage for de- 
vice B at different magnetic fields (see legend, T = 440 mK, 
V hias = 150/i V, V bg = -49.74 V, V PG and V TG 3 = 5 V). The 
inset shows as a function of magnetic field the background 
resistances that we subtracted (solid green line) and the two- 
terminal resistance with all top gates at ground (dashed 
blue line, measured at Vbias — 500 /xV). In the low field 
regime the background resistance was determined by align- 
ing the plateaus near G=4e 2 //i and in the high field regime 
by aligning the plateaus near G=8e 2 //i (b) Extracted top 
gate voltages at the plateau centres versus magnetic field. 
By fitting a linear curve to the large field data we obtain 

CVgi&2 = ^fr dv**u=A ~ L04 F/m2 - The value for Ctgi ^ 
extracted from the v — 8 data is very similar. 

dimensional ballistic channels and the formation of a 
quantum point contact. [2 lL |22| 

The 2e 2 /h steps in conductance are surprising given 
that there is both spin and valley degeneracy in bulk 
bilayer graphene, so 4e 2 /h steps are expected. The same 
observation was recently made on suspended gate-defined 
bilayer constrictions [l4[ but it remains to be understood. 
Increasing the temperature to 440 mK did not change 
the general behavior. Lowering the bias to 50 fiV (lock- 
in measurement with an AC excitation of 10 /iV) did 
not change the appearance of the plateaus either. Both 
a larger bias and a higher temperature smoothened out 
universal conductance fluctuations, as can be expected. 

We can estimate the width of the constriction from the 
position of the plateaus in top gate voltage with respect 
to the conductance minimum (Fig. [2t). We see that as 
subbands become occupied, the width of the constriction 
increases from W\ ~ 90 nm to Ws ~ 120 nm for device A 
and from W\ ~ 80 to Ws ~ 160 nm for device B, which is 
characteristic of a smooth confining potential (with W n 
the constriction width for the n-th subband). The litho- 
graphically defined separation between the respective top 
gates was less than 30 nm for both devices. This implies 
that the channel extends below the top gates, which can 
be expected given the modest band gap induced under- 
neath. 

In Fig. [3^ we explore the influence of a perpendicular 
magnetic field on the ID channels (device B). In the low- 
field regime, the plateau at 6e 2 /h quickly disappears, but 
the plateaus at 4e 2 /h and 2e 2 /h remain visible. At large 
fields, plateaus in conductance develop at 4e 2 jh (y = 4) 
and 8e 2 /h (y = 8), typical values for the quantum Hall 
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FIG. 4. (a) Current versus the voltage on TG2 (device A, 
Vtgi = Vtg3 = 10.1546 V, Vbg = -49.5 V and V hias = 
50 //V). (b) Differential conductance plotted in color scale 
versus gate voltage and source-drain bias (device A, Vtgi = 
Vtgs = 10.1546 V, Vbg = -49.5 V). The red dashed lines 
indicate the outline of a diamond as used for the analysis 
of the capacitances. The tiny diamond in the middle is not 
included in the analysis. 

effect in bilayer graphene. [23] This transition from size 
confinement to magnetic confinement occurs when the 
cyclotron radius (r c ) is equal to or larger than W n /2.[24] 
By extrapolating the positions of the plateaus for n = 2 
and v = 4, we can determine a cross-over magnetic field 
of 0.9 ± 0.2 T and estimate the size of the constriction 
based on r c (Fig. Eb). This gives W2 = 76 ± 18 nm. 
The agreement with the estimate based on the plateau 
positions (~ 120 nm) is better than a factor two. 

When we induce two barriers by appropriate gate volt- 
ages, the device behaviour changes drastically. Fig. |4^ 
shows a gate voltage scan of TG2 while TGI and TG3 are 
also biased (device A). We see sharp conductance peaks 
separated by regions of strongly suppressed conductance, 
which is characteristic of Coulomb blockade. [25] As ex- 
pected for Coulomb peaks, their position on one gate 
voltage axis varies smoothly (linearly) when another gate 
voltage is swept (Fig. [5fct) . We note that the resistance 
in Coulomb blockade is orders of magnitude larger than 
the sum of the two barrier resistances, which saturates 
around h/e 2 (Fig. [2|), indicating each gate couples to 
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FIG. 5. (a) Current plotted in color scale versus the volt- 
age applied to PG and TG3 (device B, V T gi&2 = 10.2548 
V, Vbg = -49.38 V and V hias = 50 /xV ). Occasionally 
shifts occur such as the one at Vpg ~ 10.4835 V, which 
we attribute to a rearrangement of background charges in 
the substrate, (b) Zoom-in on a conductance peak plotted 
on a logarithmic scale along with fits by a Lorentzian and 
hyperbolic cosine function (device B, Vtgi&2 = 10.1907 V, 
Vpg = 10.3880 V, Vbg = -49.86 V and V bias = 50 /xV). Fit 



^TG, off set > 



where a ■ 



2E C 



formulas are G oc cosh 2 (a- 2 _ ; , T } w .xw ,« ~~ Ar 
(the factor 2.5 assumes temperature exceeds level spacing) 

and G oc t-ft? t? — wm^ttft^, where T is the cou- 

pling strength to the leads. 



both barriers. Coulomb blockade is confirmed further by 
the diamond-shaped regions of suppressed conductance 
seen in a color plot of conductance versus gate voltage 
and bias voltage (Fig. Ht). 

The Coulomb peak spacing exhibits a standard devia- 
tion of 1.4 mV, compared to an average peak spacing of 4 
mV (device B). Irregular peak spacings can, in princi ple, 
originate from quantized level spacing contributions. [25j 
However, no excited state features are visible in the 
Coulomb diamonds so it is not clear that level spacing 
is responsible for the peak spacing variation. Instead the 
reason may be the breakdown of the constant-interaction 
model, for instance due to disorder-induced variations in 
the island size or position as charge carriers are added 
one by one, or a second island capacitively coupled to 
the main island. 

From the Coulomb diamonds and peak spacing we can 
obtain information on the dimensions and location of the 
island. For the dataset of Fig. [4)3 (device A), the addition 
energy E a dd is 0.35 ± 0.02 meV. When we neglect orbital 
energies, E a dd = 2i^ c , where E c is the charging energy. 

2 

Using E c = we can calculate the total capacitance 
of the dot, C E S = 0.46 ± 0.03 £F. From the slope of the 
Coulomb diamond edges, we find that is dominated 
by the source and drain capacitances (C s = 0.10 ± 0.024 
£F, C d = 0.30 ± 0.033 £F), which makes it difficult to 
estimate the island size from E a dd and Cs- 

Instead, we compare the measured top gate capaci- 
tances with top gate capacitances simulated using a 3D 
Poisson equation solver (Ansoft Maxwell). From the av- 
erage Coulomb peak spacing AVtg2 we extract a capac- 
itance Ctg2 ~ 45 aF. Ctg3, Ctgi and Cpg are com- 



parable. The calculated capacitance between a circular 
island of 320 nm in diameter and a metal plate with a 
320 nm hole, 50 nm above the island, is 40 aF, about as 
large as each measured top gate capacitance by itself. We 
thus infer that the island is formed in the central region 
uncovered by the top gates and extends underneath all 
the gates. This is consistent with the quantum point con- 
tacts extending underneath the split gates as discussed 
above (see also |14j). Presumably lateral confinement is 
less tight in these devices than in GaAs split gate devices, 
due to the much smaller band gap. 

Finally, we investigate whether the dot is weakly or 
strongly tunnel coupled to the leads by inspecting the 
line shape of the Coulomb peaks at low bias voltage. In 
Fig. [5Jd, we overlay on the data a fit by a hyperbolic 
cosine function and by a Lorentzian function. The for- 
mer fits the data much better, indicating that the dot 
is in the weakly coupled regime, where temperature de- 
termines the linewidth rather than tunnel coupling to 
the reservoirs. [26] Doing the same analysis for 50 peaks 
we find that they are consistently in the weakly coupled 
regime. We extract an electron temperature of 69 ± 14 
mK. 

Concluding, we show evidence for one-dimensional and 
zero-dimensional confinement in substrate supported bi- 
layer graphene split gate devices. Transport through a 
single-barrier device shows signs of conductance quanti- 
zation and in double-barrier devices clear Coulomb block- 
ade is observed. Further development of this platform 
will enable us to study the rich physics of quantum con- 
fined charge carriers in bilayer graphene. Questions to 
be addressed include how electron-electron and electron- 
phonon interactions, and spin and valley lifetimes are 
modified by the confinement. 
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